Slow-light optical bullets in arrays of nonlinear Bragg-grating waveguides 
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We demonstrate how to control independently both spatial and temporal dynamics of slow light. 
We reveal that specially designed nonlinear waveguide arrays with phase-shifted Bragg gratings 
demonstrate the frequency-independent spatial diffraction near the edge of the photonic bandgap, 
where the group velocity of light can be strongly reduced. We show in numerical simulations that 
such structures allow a great flexibility in designing and controlling dispersion characteristics, and 
open a way for efficient spatiotemporal self-trapping and the formation of slow-light optical bullets. 
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Nonlinear response is the fundamental property of op- 
tical materials that leads to interaction of propagating 
optical waves and allows all-optical switching for various 
applications. Efficiency of nonlinear interactions can be 
greatly enhanced in the regime of light propagation with 
small group velocities 1]. The slow-light regime can be 
realized when the frequency is tuned close to the edge of a 
photonic bandgap where optical waves experience strong 
Bragg scattering from a periodically-modulated dielec- 
tric structure. These effects were studied extensively in 
the structures where the propagation direction is fixed 
by the waveguide geometry, including experimental ob- 
servations of pulse propagation in fibers 2] and AlGaAs 
ridge waveguides 3] with Bragg gratings, or coupled- 
defect waveguides in photonic crystals with two- or three- 
dimensional modulation of the refractive index (see, e.g., 
Refs. 0,0,0)- I n particular, it was demonstrated that 
the group velocity and the corresponding pulse delay can 
be tuned all-optically in nonlinear photonic structures, 
and at the same time nonlinearity may compensate for 
pulse broadening due to dispersion |2|. 

All-optical control of the spatial beam dynamics be- 
comes possible in periodic two- (2D) or three-dimensional 
(3D) photonic structures, where light can propagate in 
various directions. The possibility for beam steering is 
inherently linked to the effect of beam spreading due 
to diffraction. Similar to pulses, nonlinearity can sup- 
press beam spreading and support the formation of non- 
diffracting spatial optical solitons 7]. Recent studies 
have emphasized many unique properties of spatial soli- 
tons in periodic structures such as waveguide arrays or 
photonic lattices || E| , where modulation in the trans- 
verse spatial dimension fundamentally affects the nonlin- 
ear wave dynamics. There appear multiple angular band 
gaps supporting the formation of spatial gap solitons, 
which possess tunable steering properties 10]. Recent 
experiments [llj demonstrated the effect of simultaneous 
spatial and temporal self-trapping of optical pulses in the 
form of optical light bullets |l2|. 

In this Letter, we address an outstanding key problem 
and demonstrate how to perform dynamical tunability 




FIG. 1: (a,b) Schematic of a waveguide array with (a) in- 
phase (ip = 0) and (b) phase-shifted (ip = tt) Bragg gratings; 
(c,d) corresponding isofrequency contours for different detun- 
ing from the gap edge: u) — 4 (solid), u) — 2 (dashed), uj — 1.1 
(dashed-dotted). (e,f) Corresponding diffraction patterns for 
a pulse with the central frequency at uj = 1. 



over both the magnitude and direction of the speed of 
light through all-optical control in the slow-light regime. 
In order to realize an effective nonlinear control, it is 
necessary to balance the effects of temporal dispersion in 
the propagation and spatial diffraction in the transverse 
directions. For this purpose, instead of usually studied 
symmetric photonic crystals , here we suggest to con- 
sider the photonic structures with multiple-scale modu- 
lations in different directions. This overcomes the lim- 
itations of the previously studied 2D or 3D defect-free 



photonic-crystal structures where diffraction usually in- 
creases for smaller group velocities due to a shrinkage 
of the isofrequency contours [l4|, thus greatly restrict- 
ing the efficiency of spatial self- focusing of slow light. 
Indeed, only quasi- localized nonlinear states with large 
spatial extent were predicted for a slab 2D waveguide 
with a ID Bragg grating In a sharp contrast, in this 
Letter we show that it is possible to realize the frequency- 
independent diffraction near the edge of the bandgap in 
specially designed Bragg-grating waveguide arrays. This 
allows to engineer independently the strength of diffrac- 
tion and dispersion in the slow-light regime, and thus 
provides optimal conditions for the nonlinear control of 
spatiotemporal pulse dynamics. In particular, we predict 
and demonstrate numerically the formation of strongly 
localized slow-light optical bullets in such structures. 

The arrays of nonlinear optical waveguides which are 
homogeneous in the propagation direction have been ex- 
tensively studied in recent years, and many possibilities 
for spatial beam control were demonstrated experimen- 
tally H E| . We reveal that new regimes of the spatiotem- 
poral dynamics can be achieved in the waveguide arrays 
modulated in the propagation direction with the period 
satisfying the Bragg-resonance condition at the operating 
frequency. Such structures can be fabricated in AlGaAs 
samples, with accessible fast nonlinear ity. 

In the vicinity of the Bragg resonance, the evolution of 
optical pulses can be modeled by a set of coupled-mode 
nonlinear equations for the slowly varying envelopes 
of the forward (u n ) and backward (v n ) propagating fields 
in each of n— th waveguide. In the normalized form, these 
equations can be written as follows, 



intensities) can be fully defined through the properties of 
Floquet-Bloch eigenmodes, 
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where t and z are the dimensionless time and propagation 
distance, respectively, C is the coupling coefficient for the 
modes of the neighboring waveguides, p n characterizes 
the efficiency of scattering from the Bragg grating, 7 is 
the nonlinear coefficient, and the group velocity far from 
the Bragg resonance is normalized to unity. 

We reveal that both diffraction and dispersion can be 
precisely tailored by introducing a phase shift between 
the otherwise equivalent waveguide gratings, as illus- 
trated in Figs. IHa,b). Only two- and three- waveguide 
nonlinear couplers with in-phase gratings were analyzed 
before 0, 0. Here, we consider the effect of linear 
phase shift of the gratings across the array characterized 
by the scattering coefficients p n = pexp(-ipn). With no 
loss of generality, we can take p to be real and positive. 
Then, the wave propagation in the linear regime (at small 



u n — uo exp [iKn + i/3z — tout] , 

v n = vq exp [i(K + (p)n + i(3z — iujt] , 



(2) 



where K and (3 are the transverse and longitudinal com- 
ponents of the Bloch wavevector. After substituting 
Eqs. (0) into the linearized equations (QJ (with 7 = 0), 
we obtain the following relations, 



uq[uo - (3 + 2Ccos(iT)] + v p = 0, 
u p + [w + P + 2C cos(K + if)] = 0. 



(3) 



These are the eigenmode equations which define the dis- 
persion properties of the Bloch waves, uj(K,/3). Since 
Eqs. (0) represent a square 2x2 matrix, there will ap- 
pear two branches of the dispersion curves. These de- 
pendencies determine the key features of the wave spec- 
trum. First, the gaps may appear for a certain frequency 
range, where the propagating waves with real K and 
(3 are absent. We notice that quasi-2D spectral gaps 
can indeed appear for the modes localized in the high- 
index waveguides, i.e. below the light-line of the sub- 
strate. Second, the spatial beam refraction and diffrac- 
tion are defined by the shape of the corresponding isofre- 
quency contours. The propagation angle is defined by the 
value of a(cj,K) = —d/3/dK, and the effective diffrac- 
tion experienced by the beam depends on the curvature, 
D(uj,K) = —d 2 (3/dK 2 . Below, we analyze how these 
fundamental characteristics can be controlled by select- 
ing the phase shift cp between the Bragg gratings. 

For the in-phase gratings [i.e., when <p = 0; see 
Fig. HJa)], the dispersion relation becomes uj{K,f3) — 
-2Ccos(K) ± (p 2 + (3 2 ) 1 ' 2 . It follows that the 2D 
gap appears only when the waveguide coupling is weak, 
C < p/2. On the other hand, the shape of the isofre- 
quency contours is strongly frequency-dependent as the 
transmission band edge is approached, see Fig.^c). This 
happens because the position of the one-dimensional fre- 
quency gap depends on the propagation direction. 

We find that the nature of wave dispersion is funda- 
mentally altered for the waveguide structure with out- of - 
phase shift of the neighboring gratings [i.e., when tp = 7r; 
see Fig.QJb)]. The corresponding dispersion relation has 
a different form, lj(K,0) = ±{p 2 + [/? - 2Ccos(iT)] 2 } 1 / 2 . 
Moreover, for any propagation angle defined by the trans- 
verse Bloch wavevector component the width and po- 
sition of the one-dimensional frequency gap remains the 
same, \uj\ < p. This unusual property leads to remark- 
able spectral features. First, the 2D (quasi-) gap is al- 
ways present in the spectrum irrespectively to the grat- 
ing strength (p) and coupling between the waveguides 
(C). Second, the shape of isofrequency contours does 
not depend on frequency in the transmission band, see 




FIG. 2: (colour online) Family of the stationary (zero- 
velocity) light bullets characterized by (a) energy and 
(b) width along the transverse (solid) and longitudinal 
(dashed) directions vs. the frequency tuning inside the spec- 
tral gap. (c-d) Intensity profiles of slow- light optical bullets 
for (c) u) — 0.8 and (d) u) — 0.995; brighter shading marks 
higher intensity. 
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FIG. 3: (colour online) Moving slow- light optical bullets 
with the frequency detuning co = 0.5. Notations are the same 
as in Fig. and the velocity is normalized to the speed of 
light away from the Bragg resonance. Profiles are shown for 
the velocities (c) V = 0.1 and (d) V = 0.2. 



Fig.njd). This means that the beam diffraction remains 
the same even when the band edge is approached. 

The dependence of diffraction on frequency strongly 
affects the shaping of short pulses with broad frequency 
spectra. We perform numerical modeling of the pulse 
dynamics, when the input beam is focused into a single 
waveguide. The pulse duration is T = 10, and the central 
frequency is ujq = 1. The snapshots of the light intensity 
after the pulse enters the photonic lattice are shown in 
Figs.EJe,f). For the structure with the in-phase gratings 
[Figs, me)], the spatial diffraction strongly depends on 
frequency, creating a colored pattern similar to the su- 
perprism effect in photonic crystals [HI- For the phase- 
shifted gratings, all the components inside the frequency 
band experience exactly the same diffraction and propa- 
gate together. As a result, the intensity profile [Figs-GJe)] 
has precisely the same shape as for the discrete diffrac- 
tion of monochromatic beams 8], with well-pronounced 
field minima and peak intensities at the outer wings. 

The unique features of linear spectrum in arrays with 
phase-shifted gratings suggest that these structures pro- 
vide optimal conditions for a nonlinear control of the 
pulse dynamics. In particular, since the 2D gap ap- 
pears for any values of the grating strength and waveg- 
uide coupling, it is possible to choose these parameters 
independently in order to balance the rates of dispersion 
and diffraction. This allows for simultaneous compensa- 
tion of the pulse broadening in space and time through 
the nonlinear self-trapping effect. Indeed, we find nu- 
merically localized solutions of Eqs. for stationary 
and moving light bullets of the form of solitary waves 
{u, v} n (z, t) = {2, v} n (z — Vt) exp(— iwi), where V is the 
propagation velocity. We confirm that localization is pos- 
sible in both the cases of positive (7 = +1) and negative 
(7 = —1) nonlinear response, since anomalous or nor- 
mal dispersion regimes can be accessed on either edges 
of the photonic bandgap. The soliton solutions for dif- 
ferent signs of 7 can be mapped by changing uo —> —uo 
and making a corresponding transformation of pulse pro- 
files. In Fig. [21 we show the characteristics of immobile 
solitons such as energy and width vs. the frequency de- 
tuning inside the bandgap for self-focusing nonlinearity. 
These solitons have a well-pronounced X-shape near the 
upper edge of the gap, becoming more localized inside 
the gap as the pulse energy is increased. We notice how- 
ever that these are fully localized states that should not 
be confused with the so-called X-waves |2£| which re- 
main quasi-localized only over a finite propagation dis- 
tance. The gap solitons can also propagate along the 
waveguides, and we present the characteristic of moving 
solitons in Fig. These slow-light bullets become more 
extended in both the transverse and longitudinal direc- 
tions as the propagation velocity is increased. 

Finally, we perform numerical simulations of the spa- 
tiotemporal pulse dynamics in these structures. In the 
linear regime, the pulse broadens in both transverse and 
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FIG. 4: Snapshots of field intensities for an optical pulse 
propagating in a waveguide array structure show in Fig. Qd) : 
(a-d) Linear broadening due to spatial diffraction and tempo- 
ral dispersion; (e-h) Nonlinear self-trapping in space and time 
and formation of an optical bullet. 



longitudinal directions [Figs. Sfa-d)]. Nonlinear self- 
action results in the pulse self-trapping in both space 
and time [Figs. 0f e-h)]. In this example, the velocity of 
the generated light bullet is below 30% of the speed of 
light in the absence of the Bragg grating, and smaller 
velocities can be accessed as well by controlling the fre- 
quency of the input pulse. The profile of the self-trapped 
state shown in Fig. Sfh) closely resembles the shape of 
the exact soliton solution shown in Fig. Eld). 

In conclusion, we have revealed that both spatial 
diffraction and temporal dispersion can be engineered in- 
dependently near the edge of the photonic bandgap in 
the waveguide arrays with phase-shifted Bragg gratings. 
We have shown that such structures possess quasi-two- 
dimensional photonic bandgaps and provide optimal con- 
ditions for self-localization of pulses in nonlinear media. 



We have demonstrated that these waveguide arrays can 
be employed for shaping and control of optical pulses si- 
multaneously in space and time, and allow the formation 
of optical bullets propagating with slow group velocities. 
Such slow-light optical bullets offer novel possibilities for 
all-optical switching, steering, and control of short pulses. 
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tralian Research Council. 
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